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ABSTRACT 

We present X-ray analysis of two Wolf-Rayct (WR) binaries: V444 Cyg and CD 
Cru using the data from observations with XMM-Newton. The X-ray light curves show 
the phase locked variability in both binaries, where the flux increased by a factor of ~ 2 
in the case of V444 Cyg and ~ 1.5 in the case of CD Cru from minimum to maximum. 
The maximum luminosities in the 0.3-7.5 keV energy band were found to be 5.8 x 10 32 
and 2.8 x 10 32 erg s _1 for V444 Cyg and CD Cru, respectively. X-ray spectra of these 
stars confirmed large extinction and revealed hot plasma with prominent emission line 
features of highly ionized Ne, Mg, Si, S, Ar, Ca and Fe, and are found to be consistent 
with a two-temperature plasma model. The cooler plasma at a temperature of ~ 0.6 
keV was found to be constant at all phases of both binaries, and could be due to a 
distribution of small-scale shocks in radiation-driven outflows. The hot components 
in these binaries were found to be phase dependent. They varied from 1.85 to 9.61 
keV for V444 Cyg and from 1.63 to 4.27 keV for CD Cru. The absorption of the hard 
component varied with orbital phase and found to be maximum during primary eclipse 
of V444 Cyg. The high plasma temperature and variability with orbital phase suggest 
that the hard-component emission is caused by a colliding wind shock between the 
binary components. 

Key words: stars:Wolf-Rayet - stars:binary - stars:massive - stars:individual(V444 
Cyg, CD Cru) - starsX-rays 



1 INTRODUCTION Trapezium cluster (Stelzer et al. 2005), M17 cluster (Broos 

et al. 2007) and Carina Nebula (Leutenegger & Kahn 2003) 
WR stars are known to produce strong stellar winds driven . jj-j.-it.jv i • • 

° require an additional hard X-ray component to explain their 

by their strong radiation held. The stellar winds can reach „ . , , , , , . . 

J ° _, X-ray spectra, this additional hot component is mconsis- 

velocities up to 1000-3000 km s with the mass loss rates .., n , , , n . n , n ~ 

_ fi _, tent with the standard model which produces only a sort 

of 10 — 10 Mo yr depanding upon mass and age v T , n . ,, , , , „ 

. X-rays. In order to explain the hard X-ray component, Ba- 

(Rauw 2008). These winds not only afreet their evolution, , , „ ,, , ,.„„_> , - , ,. „ n , 

v ' J bel & Montmerle (1997) put forward a magnetically confined 

but also have a tremendous impact on their surroundings , , , n n . ,. ... , n . „ . , 

° wmd-shock model. According to this model, the stellar wind 

(Rauw 2008). Most of the O-type and early B-type single . „ j , n , ,. „ ,j . . 

v ' i, o„ _i N . is conhned by the large-scale magnetic held into the equa- 

stars arc reasonably bright (10 < Lx < 10 erg s ) and ... ... j i 

J ° y , tonal region and the two streams from the upper and fower 

are soft (kT < 1 keV) X-ray sources. The X-ray luminos- , . n ,,. , ., n ... , , , . , 

v ' J J hemispheres collide with it to heat the plasma to temper- 

ity of these sources arc found to scale with the bolomet- , , . , . „ , , r . n 

J . _ 7 atures much higher than typically expected tor singfe stars 

ric luminosity as Lx/Lboi ~ 10 (Berghofer et af. 1997, 



Sana et al. 2006) and is generally thought to be produced 



from the standard wind-shock model. However, the magnet- 
ically confined wind-shock model explains hard X-ray on 



due to the shocks, with velocity iumps up to a few hun- , „ , r . , , , n \ r 

only a small number ot single early-type stars, e.g.,6 1 On C 



(Gagne et al. 2005) 



dred km s , generated throughout the stellar wind due to 
dynamic instabilities (Lucy & White 1980; Owocki & Co- 
hen 1999; Kudritzki & Puis 2000). This 'wind-shock' see- Stellar binary systems with early type (WR or OB) 
nario has now become the 'standard' model to explain the components are found to be more luminous in X-ray than 
X-ray emission from single early-type stars (Rauw 2008 and what is expected from the individual components separately 
references therein). But several massive stars in the Orion (De Becker 2007). This excess X-ray luminosity is thought to 
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Table 1. Stellar parameters of WR binaries V444 Cyg and CD Cru. 



Parameters 


V444 Cyg 


Ref a 


CD Cru 


Ref a 


Other Names 


HD 193576; WR 139 




HD 311884 ; WR 47 




Spectral Type 


OR TTT-V -4- WN5 


1 


05 V 4- WN6 


1 


Mass (M Q ) 


25 + 10 


2 


57 + 48 


10 


V (mag) 


7.94 


3 


10.89 


11 


Distance (pc) 


1900 


3 


3000 


12 


A v (mag) 


2.48 


3 


3.56 


12 


Orbital Period (days) 


4.21 


2 


6.24 


13 


Binary separation (R0) 


38 


2 


68 


10 


inclination (°) 


78 


4 


70 


10 


eccentricity 


0.03 


5 





11 


M (lCT 5 M Q yr" 1 ) 


0.06 + 0.60 


6,4 


0.10 + 3.00 


15,10 


Voo (km s _1 ) 


2540 + 1785 


7.8 


3000 + 2460 


15,8 


log(Lboi) erg s" 1 


39.24 + 38.56 


9 


39.53 + 38.69 


9 


Radius (R©) 


10 + 3 


5 


57 + 48 


10 


Cluster Membership 


Berkeley 86 


3 


Hogg 15 


11 



References: 1. van der Hucht et al. (1981) 2. M finch (1950) 3. Massey, Johnson & Dc Gioio-Eastwood (1995) 4. Kurosawa., Hillier & 
Pittard (2002) 5. Cherepashchuk, Eaton & Khaliullin (1984) 6. Marchenko et al. (1994) 7. Underhill & Fahey (1987) 8. Nugis & Lamers 
(2000) 9. Landolt-Bornstein (1982) 10. Moffat et al. (1990) 11. Moffat (1974) 12. Sagar, Munari & Boer (2001) 13. Nicmela, Massey, 
Conti (1980) 14. Cherepashchuk & Karetnikov (2003) 15. Rcpolust, Puis & Herrero (2004) (typical value for the spectral-type) 



originate from the wind-wind interaction zone where post- 
shock temperature is believed to reach up to a few to sev- 
eral tens of million Kelvin (MK), contributing mainly to 
the harder energy (kT > 1 keV) X-ray emissions. The post- 
shock temperature is expected to be higher in long period 
(~ a few weeks) binaries as the stellar winds collide close 
to their terminal velocities. In short period (~ a few days) 
binaries, the wind collision occurs while winds are still accel- 
erating (De Becker 2007). Though, in both the cases, models 
predict phase-locked variations in the X-ray domain (Luo, 
McCray & Mac Low 1990; Usov 1992; Antokhin, Owocki 
& Brown 2004). These variations are produced either due 
to the varying circumstellar optical depth along the line of 
sight towards the shock as the stars revolve around each 
other or by varying the orbital separation in eccentric bi- 
naries, which changes the intrinsic strength of the collision 
(Rauw 2008; Naze et al. 2007). On the contrary, there are ex- 
amples of a number of massive star binaries (O+O) which 
show a weak enhancement in the X-ray emission over the 
single massive stars, but present no variability in the X-ray 
flux (Albacete Colombo, Mendez & Morrell, 2003). Antokhin 
(2007) showed that, in general, the X-ray properties of single 
and binary early- type stars are not very different. 

In this paper, we investigated the X-ray emission char- 
acteristics for a two short-period WR binaries, namely V444 
Cyg and CD Cru using high signal-to-noise (S/N) ratio 
XMM-Newton observations. For both objects, the X-ray 
data from XMM-Newton is analysed for the first time, while 
for CD Cru, we analyze the X-ray emission for the first time. 
Both the stars are collected from a systematic study car- 
ried out on young galactic clusters using the XMM-Newton 
archival data. The choice of these stars associated with 
young clusters provide us a good knowledge of the distance 
and the reddening parameter. Details of the complete study 
will be presented elsewhere. The relevant parameters of the 
sample are listed in Table [T] and a description on the indi- 
vidual objects are given below. 

V444 Cyg (WR 139) is an 06 III-V + WN5 eclips- 



ing binary system associated with the cluster Berkeley 86 
(van der Hucht 2001) situated at an estimated distance of 
1.9 kpc with E(B-V)~0.8 mag (Massey, Johnson & Degioia- 
Eastwood 1995). The orbital period and an eccentricity of 
the binary system are estimated to be 4.21 days (d) and 
0.03, respectively, (Khaliullin, Khaliullina & Cherepashchuk 
1984). It has been observed in multi- wavelengths with well 
determined physical parameters. A fairly soft (kT = 0.5 keV) 
thermal X-ray spectrum was reported using Einstein obser- 
vations (Moffat et al. 1982). The flux was found to vary by 
a factor of 2 with a minimum at phase zero, which implies 
that the WN5 star is in front of the 06 star. Pollock (1987) 
reported an absorption-corrected flux of Lx = 7.7 x 10 32 erg 
s _1 in the 0.4-4 keV range. In ROSAT observations, Corco- 
ran et al. (1996) reported a minimum absorption-corrected 
X-ray luminosity of 3.8 x 10 32 erg s _1 at the primary eclipse 
and 8.6 x 10 32 erg s _1 at the secondary eclipse in the 0.4- 
4.5 keV range. The ASCA observations at three phases in 
the energy band of 0.7-10 keV were analyzed by Maeda et 
al. (1999). They found the presence of two-temperature (0.6 
and 2.0 keV) plasma, with soft component (0.2-4 keV) lu- 
minosity of 6 — 11 x 10 32 erg s _1 . The high temperature 
plasma (kT ~ 2 keV) and the variability with the orbital 
phase suggest that the hard component emission is caused 
due to the shocks originated from the collision of winds from 
WN5 and 06 stars (Maeda et al. 1999). 

CD Cru (WR 47) is an 05 V + WN6 binary system 
reported as a member of the cluster Hogg 15 (van der Hucht 
2001). Throughout this study, we adopted a distance of 3 
kpc of Hogg 15 (see Piatti et al. 2002 ;Sagar, Munari & 
Boer 2001). The orbital period of this binary system is 6.24 
d with zero eccentricity (Niemela, Massey & Conti 1980). 
No previous X-ray observations of the system are reported 
in the literature. 

Our paper is organised in the following manner; in Sec- 
tion 2 and 3, we give a detailed description of X-ray data 
reduction and analysis, respectively. In Section 4, we present 
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Table 2. Journal of XMM-Newton Observations of the objects in the sample. 



Object Name 


Observation 


Exposure 


Start time 


EPIC filter 


Offset from 


Orbital 




ID 


Time (sec) 


UT(hh:mm:ss) 




Target (arcmin) 


Phase 


V444 Cyg 


0206240201 


10032 


19 May 2004 12:35:42 


Thick 


0.001 


0.13 (0.12-0.14) 




0206240301 


11914 


27 May 2004 12:01:39 


Thick 


1.120 


0.03 (0.01-0.04) 




0206240401 


11915 


29 May 2004 12:01:49 


Thick 


1.162 


0.51 (0.50-0.51) 




0206240501 


11918 


06 Juno 2004 11:36:30 


Thick 


1.111 


0.39 (0.38-0.41) 




0206240701 


11911 


14 June 2004 10:58:29 


Thick 


1.133 


0.29 (0.27-0.30) 




0206240801 


19921 


27 Oct 2004 23:47:44 


Thick 


1.093 


0.47 (0.45-0.50) 


CD Cru 


0109480101 


53046 


03 July 2002 15:51:56 


Thick 


2.812 


0.78 (0.73-0.83) 




0109480201 


48015 


26 Aug 2002 21:55:14 


Thick 


3.848 


0.47 (0.43-0.51) 




0109480401 


53040 


21 Jan 2003 01:07:38 


Thick 


4.969 


0.05 (0.01-0.09) 



a discussion on the X-ray properties of the sources. Finally, 
we summarize and draw our conclusions in Section 5. 



2 X-RAY OBSERVATIONS AND DATA 
REDUCTION 

The log of the X-ray observations is given in Table [2] We 
made use of the archival data obtained with the XMM- 
Newton observatory which consists of three co-aligned X-ray 
telescopes observed simultaneously and covered 30' x 30' re- 
gion of the sky. The X-ray photons were recorded with the 
European Photon Imaging Camera (EPIC), which forms im- 
ages on three CCD-based detectors: the PN (Stritder et al. 
2001), and the twin MOS1 and MOS2 (Turner et al. 2001) 
with an angular resolution of 6" (FWHM). During the ob- 
servations, all the three EPIC detectors were active in full 
frame mode together with the Thick filter. For V444 Cyg, six 
separate observations were taken with exposure time rang- 
ing from 10 to 20 ks. These observations are spread over the 
half cycle of binary system. However, for the star CD Cru, 
three separate observations were taken covering almost full 
binary cycle. 



2.1 EPIC data reduction 

We reduced the X-ray data using standard XMM-Newton 
Science Analysis System software (SAS version 7.0.0) with 
updated calibration files (Ehle et al. 2004). Event files for 
the MOS and the PN detectors were generated using the 
tasks emchain and epchain, respectively. These tasks allow 
calibration of the energy and the astrometry of the events 
registered in each CCD chip and to combine them in a sin- 
gle data file. We restricted our analysis to the energy band 
0.3-7.5 keV, as the data below 0.3 keV are mostly unrelated 
to bonafide X-rays, and above 7.5 keV is mostly dominated 
by the background counts. Event list files were extracted us- 
ing the SAS task evselect. Data from the three cameras 
were individually screened for the time intervals with high 
background when the total count rates (for single events of 
energy above 10 keV) in the instruments exceed 0.35 and 
1.0 counts s _1 for the MOS and PN detectors, respectively. 
The observations with observation ID 0206240401 for the 
source V444 Cyg were heavily affected by the high back- 
ground events, and the data in MOS detectors were not use- 
ful. Only half of the observation time (~5 ks) of PN was 
found useful. 



Light curves and spectra were extracted using a circu- 
lar region with the source as the center in the energy range 
0.3-7.5 keV of the EPIC detectors. The X-ray sources in 
the cluster were often found to be largely contaminated due 
to the emission from the neighboring sources. For this rea- 
son, the radii of extraction regions were varied between 20" 
and 30" depending upon the position of the sources in the 
detector and their angular separation between the neighbor- 
ing X-ray sources. The background has been estimated from 
a number of empty regions close to the X-ray source in the 
same CCD of the detector. X-ray spectra of the sources were 
generated using SAS task especget, which also computes 
the photon redistribution as well as the ancillary matrix. For 
each source, the background spectrum was obtained from re- 
gions devoid of any sources chosen according to the source 
location. Finally, the spectra were re-binned to have at least 
20 counts per spectral bin for both the sources in all the 
observations. 

2.2 RGS data reduction 

The Reflection Grating Spectrometers (RGS) are mounted 
on two of the three XMM-Newton X-ray telescopes and 
were operated in spectroscopy mode during the observa- 
tions. We followed the standard procedure as outlined in 
the XMM-Newton handbook to generate the RGS spectra 
of the sources. The raw data were processed with the task 
RGSPROC at the position of the sources. The other bright 
sources, found in the same field observed by RGS, have been 
excluded from background estimation. In addition, we have 
filtered the event list for high background level events. For 
star V444 Cyg, the data obtained in observations with ob- 
servation ID 0206240401 were contaminated with high back- 
ground events in RGS detector, therefore the data was not 
useful. For CD Cru, the two observations ID 0109480101 
and 0109480201 were not useful as they are heavily contam- 
inated by a brighter source HD 110432 and affected by high 
background episodes, respectively. 



3 ANALYSIS AND RESULTS 
3.1 X-ray light curves 

The background subtracted X-ray light curves of the WR- 
binaries V444 Cyg and CD Cru, as observed with PN detec- 
tor are shown in Fig. HI a) andQJb), respectively. The light 
curves are in the 0.3-7.5 keV (total) energy band show the 
variability which suggests colliding wind shocks. Further, 
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Figure 1. X-ray light curves of the WR-binaries (a) V444 Cyg and (b) CD Cru in the 0.3-7.5 keV energy band. 
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Figure 2. (a) PN and (b) MOS light curves at three energy bands: the total (0.3-7.5 keV), the soft (0.3-2.0 keV) and the hard (2.0-7.5 
keV), and hardness ratio HR curve as a function of orbital phase, where HR = hard/soft of V444 Cyg. 



we performed the \ 2 test to measure the significance of the 
deviations from the mean count rate in order to quantify 
the constancy of the data over the time-scale of observa- 
tions. We found the variability in the light curves with a 
confidence level of greater than 99.999% for both binaries. 
In order to investigate the variability in the different en- 
ergy bands, the light curve of these binaries obtained with 
the MOS and PN data are generated into two energy bands 
namely the soft (0.3-2.0 keV) and the hard (2.0-7.5 keV). 
The hardness ratio (HR) is defined by the ratio of hard to 
soft band count rates. The total, hard and soft band in- 
tensity curves, and the HR curve as a function of the or- 



bital phase are shown in the subpanels running from top to 
bottom in Fig. [gji (PN) and[2p (MOS) for the WR binary 
V444 Cyg. Similar plots of intensity and HR for the star CD 
Cru are shown in Fig. [3] The phases of the observations are 
reckoned using the ephemeris HJD = 2441164.337+4.213E 
for V444 Cyg (Underbill, Grieve & Louth 1990) and HJD 
= 2443918.4000+6.2399E for CD Cru (Moffat et al. 1990, 
Niemela, Massey & Conti 1980). Here, the phase 0.0 indi- 
cates the primary eclipse and the phase 0.5 indicates the 
secondary eclipse. The light curves in the individual bands 
show the phase locked variability for both binaries. For V444 
Cyg, the count rates in total energy band were minimum at 
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Figure 3. Similar to Fig.[2]but for CD Cru. 



phase 0.0 and maximum at the phase 0.45. After the phase 
0.45, the count rates were decreased upto the phase 0.5. The 
count rates are increased by a factor of ~ 2 being minimum 
at phase 0.0 in the total energy band. In the soft band, 
the count rates were minimum at phase 0.0 and maximum 
during the phase 0.45-0.5. However, in the hard band light 
curve two minima at phase 0.0 and 0.5 were seen clearly. 
During the phases from 0.13 to 0.45 (i.e. outside the eclipse) 
the count rates in the hard band were constant. The HR 
curve can reveal the information about the spectral varia- 
tions. The HR curve shows that the emitted X-ray is harder 
just after phase 0.0 being maximum at phase 0.13, and af- 
terward decreased till the phase 0.5. In the case of CD Cru, 
the nature of the variability was found to be similar to that 
seen in V444 Cyg. In soft band, the count rates were found 
to be constant during the phases 0.47 and 0.78 (i.e. outside 
the eclipse). However, the count rates are decreased rapidly 
after the phase 0.47 in the hard band. A small variation 
was also observed in the HR curve of the CD Cru. It was 
maximum at the phase 0.0 and minimum during the phase 
0.47. 



3.2 X-ray spectra and spectral fits 

The EPIC spectra of WR binaries in different phases are 
shown in Figs. 2] and [S] for V444 Cyg and CD Cru, respec- 
tively. Below 1 keV, the spectra were found to be affected 
by the high extinction as previously observed with ASCA 
for V444 Cyg (Maeda et al. 1999). Strong emission lines are 
seen in the MOS and PN spectra of both WR binaries. The 
most prominent lines found in the spectra along with their 
laboratory energies are the following: Fe XVII (0.8 keV), Ne 
X (1.02 keV), Mg XII (1.47 keV), Si XIII (1.853 keV), S XV 
(2.45 keV), Ar XVII (3.12 keV), Ca XIX+XX (3.9 keV) and 
Fe XXV (6.63 keV). 

In order to trace the spectral parameters at different 
binary phases, we performed spectral analysis of each data 



set using simultaneous fitting of EPIC data by two models, 
(a) plane-parallel shock model (vpshock; Borkowski, Lyerly 
& Reynolds 2001), and (b) models of Astrophysical Plasma 
Emission Code (APEC; Smith et al. 2001), as implemented 
in the XSPEC version 12.3.0. Ax 2 - minimization gave the 
best fitted model to the data. The presence of interstellar 
material along the line-of-sight and the local circumstellar 
material around the stars can modify the X-ray emission 
from massive stars. We corrected for the local absorption 
in the line-of-sight to the source using the photoelectric ab- 
sorption cross sections according to Balucinska-Church & 
McCammon (1992) and modeled as phabs (photoelectric 
absorption screens) with two absorption components, N|f M 
and Njj cal . The Ng M was estimated using the relation, Nh 
= 5.0xl0 21 x E(B-V) cm" 2 (Vuong et al. 2003), where 
E(B — V) = Ay/3.1, assuming a normal interstellar redden- 
ing law towards the direction of the cluster. We used the val- 
ues of Av nearly 2.48 mag and 3.56 mag derived for the clus- 
ter Berkeley 86 (Massey, Johnson & Degioia-Eastwood 1995) 
and Hogg 15 (Sagar, Munari & Boer 2001), respectively. The 
estimated values of Njf M towards V444 Cyg and CD Cru 
are found to be 4.0 x 10 21 cm" 2 and 6.0 x 10 21 cm -2 , re- 
spectively. The Nn cal was estimated by making a fit to the 
observed spectra by varying the local environment for the 
soft (kTi) and the hard (kT2) energy components in terms 
of Njj and N H , respectively. Because the WN stars are at 
evolved stages, the abundances of He, C, N, O, Ne, Mg, Si, 
S, Ar, Ca, and Fe were allowed to vary during the fitting 
procedure to account for the observed line emission. The 
solar abundances were adopted from Lodders (2003). 

First, we fitted VPSHOCK plasma model to derive 
their spectral features. The constant temperature VPSHOCK 
plasma model was considered without incorporating mass- 
loss and orbital parameters of WR stars. However, the model 
does account for non-equilibrium ionization effects and as- 
sumes an equal electron and ion temperature. The best fit 
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Figure 4. X-ray spectra of MOS and PN data with the best fit 2T VAPEC model in upper subpanels of each graph for V444 Cyg. The x 2 
distribution in terms of ratio are given in lower subpanels of each graph. The last three digits of observation ID and the corresponding 
phases of the observations are given in each graph. 



VPSHOCK models to the data for the WR stars are given in 
Tableland Table|4]for V444 Cyg and CD Cru, respectively. 

X-ray emitting plasma may not be isothermal and the 
observed X-ray spectrum may be a superposition of a cool 
stellar component and a hot colliding wind shock plasma, as 
a number of emission lines in the spectra can be formed over 
a range of temperatures. The cooler component is believed 
to arise from the instabilities in radiation-driven outflows. 



However, using the X-ray imaging data of XMM-Newton 
it is not possible to resolve the X-ray emission from col- 
liding wind region and individual stars separately. There- 
fore, secondly, we fitted two temperature (2T) plasma model 
"vapec" to characterize such components. The form of 2T 
plasma model was phabs(phabs*VAPec+phabs*vapec). In 
terms of x 2 the 2T plasma model provides better goodness- 
of-fit for both the sources over the VPSHOCK model. The 
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Figure 5. Same as Figure l4l but for CD Cru. 



results for the 2T plasma model along with the data are 
displayed in Fig. g] and Fig. [S] for V444 Cyg and CD Cru, 
respectively. The best-fit parameters are given in Table [3] 
and Table Q] for V444 Cyg and CD Cru, respectively. 



Table 3. Spectral parameters derived from the best fit VPSHOCK model and two-temperature VAPEC plasma models for V444 Cyg. 



00 



V444 Cyg 


Obs ID 


0206240201 


0206240301 


0206240401 


0206240501 




0206240701 


0206240801 


Phase 


(0.12-0.14) 


(0.01-0.04) 


(0.50-0.51) 


(0.38-0.41) 




(0.27-0.30) 


(0.45-0.50) 


Model 


vpshock 2Tvapcc 


vpshock 2Tvapcc 


vpshock 2Tvapcc 


vpshock 2Tvapcc 


vpshock 


2Tvapec 


vpshock 2Tvapec 



kT 2 

EM! 

EM 2 

He 

C 

N 

O 

Ne 

Mg 

Si 



0.23 
0.19 

763.00 

9998.91 

650.45 

191.35 

68.34 

0.00 

0.00 

999.97 

6.02 

9991.50 

92.50 

2.77(130) 



0.48 
11.86 
0.57 



+ 0.20 
0.12 
+ 1.66 



+ 0.07 
0.13 
„_ + 0.22 
-0.26 



5 10+ ' 72 

< 30.15 

< 693.28 

< 651.73 
377.3S_ 33 gg 

< 80.79 
„+130.88 



26.32 
+ 152.45 

40.66 
+ 135.50 
17.19 



31.77 
62.63 
38.33 

< 119.86 

< 11.46 
20.95+ 2 
1.07 (127) 
32.691 
31.739 
32.639 



0.18 
0.17 

518.90 

0.00 

0.00 

2.13 

33.83 

76.34 

142.88 

199.06 

6.02 

2994.71 

196.91 

2.38(79) 



64+ 016 
u - o4 -0.24 

9 35 + 4 ' 23 
y ' J -2.50 

60+ 04 
UDU -0.04 
R + 5.35 
-4.16 



8.16^ 



< 0.14 
S3 0.00 
S3 0.00 

< 66.17 

< 20.27 
+ 365.20 



38.42 
37.36 
73.74 

< 64.00 

< 97.37 
66.99+!^ 
1.28 (77) 
32.485 
31.746 
32.398 



7.61 
+ 344.36 

9.20 
+ 664.80 
19.03 



+ 1.46 
'-0.18 



, + 0.74 
"-0.48 



0.38 
< 4.25 



+ 0.82 
0.19 



< 9997.63 
S3 0.00 

< 9998.63 

< 25.84 

< 11.42 
+ 289.49 

17.30 
+ 454.70 
35.98 
+460.11 
133.78 
SB 0.00 

< 995.93 

, , ,-+113.12 
14 - 27 _9.40 
1.10 (27) 



40.46 
78.56 
217.19 



0.47 



+ 0.26 
0.68 



<12.20 

62 + 012 
u o -0.13 

1 97+ 5 ' 70 
1 - -0.95 

<288.07 
<36.44 
S3 0.00 
S3 0.00 
S3 0.00 
S3 0.00 
< 150.67 

+ 225.47 



10.36 
+ 220.92 

13.96 
+ 491.05 

19.90 



29.21 
97.03 
58.63_ 
> 2.77 

< 337.62 

„ 22 +157.32 

-8.12 
0.77 (25) 
32.589 
32.170 
32.380 



„ + 0.30 
D — 0.18 



3 14+ - 92 
31 -0.43 

Q .3+O.23 
u " ' -0.20 
. 7 ,+0.65 
' -0.80 



< 951.62 
S3 0.00 

S3 0.00 

1114 90+ 966715 

< 58.76 
„ + 9.53 

-4.6B 
21. 20"^ 



5.34 



25.91_ 
S3 0.00 
89.45 



+ 56.08 

3.95 
+ 85.94 
-8.20 



+ 125.95 
-70.88 
,, ,.,+105.67 
^■ di -13.01 
1.47 (200) 



o ss + 0-20 
°' S5 -0.25 

2 74+ 1 ' 71 
2 - 7 -0.91 
„ _„ + 0.09 
°- 5S -0.10 
3.64 



+ 1.01 
0.97 



29.331^ 

< 74.03 

< 728.99 

< 258.19 
S3 0.00 

< 191.51 
10.13 
2.83 
3.19 

< 72.61 
28 



+ 64.95 
3.96 
+ 15.40 

2.13 
+ 80. SO 
5.29 



1.93"' 



+ 32.02 
-16.22 
< 37.60 
, + 12.98 
-1.45 
1.38 (198) 
32.764 
32.104 
32.657 



,,+0.15 
■ ' -0.40 



O + 0.55 
-0.46 



0.58 
< 35.92 



+ 0.18 
0.14 



617 61 +9380.65 
-201.57 

S3 0.00 

> 400.059 

3054 69 + 7 682.71 

3uo4Da — 48.20 

< 61.40 

< 46.27 
73.63j 
212. 54^ 
S3 0.00 

«n ,-+9999,82 
631 16 -260.60 
1R , ..+9836.23 

1.13 (168) 



, + 37.30 
-24.93 
,, + 109. 78 
-33.55 



8.40 
< 0.17 
0.70 
9.61 



+ 3.42 
2.88 



+ 0.15 
0.10 
+ 3.51 
1.97 



, .,+20.73 
^■ a± -2.52 

< 25.367 
> 40.18 
S3 0.00 

S3 0.00 

< 316.85 

< 446-78 

< 658-42 
120.13 
203.53 



+411.08 

68.09 
+ 250.52 
91.46 



246.82 + 380-31 
> 182.96 

142 33+ 124 ' 11 
14 2 .33_ 57 97 

1.08 (166) 

32.726 

31.999 

32.636 



, + 0.15 
-0.19 



,+0.70 
-0.78 



0.35 
< 11.77 



+ 0.10 
0.07 



QQ73 b + 9520.95 
-199.229 

< 9999.97 

< 98.928 

< 37.34 

< 2.96 
,+6.66 



4.49 
+ 16.25 
9.83 

+ 855.09 
40.00 



13.17 

47.98 

184.70 

S3 0.0 

< 9997.92 

11 35 + " 76 
11 do -7 .19 

1.12(277) 



n 78+ 012 
°- 78 _0.11 
, ,,+0.76 
^• JJ _0.84 

60 + 03 
-0.03 

3.22 



+ 0.68 
0.34 



30.32 2 



, + 50.70 
-22.53 

> 20.04 

< 109.71 

< 885.30 

< 18.09 

< 0.49 
+ 2.98 



1.43 
4.59 
44.18 
23.98 
4.89 
1.11 



0-62 
+ 33.75 
1.00 
+ 38.89 

3.08 
+ 67.11 
3.75 
+ 35.39 

4.28 
+ 4.37 
0.31 
0.98 (275) 
32.732 
32.177 
32.591 



to 



Note: Here x£ = 
and are 
kT x and kT 2 an 
t is in the units 
EMi and BM 2 a 
Lj^ luminosity ii 



X /v, where u is deg: 
in the units of 10 22 c 
; in the units of kcV. 



i of freedom. Errors arc with 90% < 



■e in the units of 10 04 cm~'\ 
band 0.3—7.5 kcV. luminosity : 



band 0.3-2.0 kcV and L" 1 



band 2.0-7.5 kcV i 



i the units of ergs 



XMM-Newton Observations of WR binaries 9 



Table 4. Spectral parameters derived from the best fit VPSHOCK model and two-temperature VAPEC plasma models for CD Cru. Units 
of the parameters are same as given in Table [3] 



CD Cru 



Obs ID 0109480101 0109480201 0109480401 

Phase (0.73-0.83) (0.43-0.51) (0.01-0.09) 



Model 


vpshock 


2Tvapec 


vpshock 


2Tvapec 


vpshock 


2Tvapec 


N H 


. ..+0.14 
1 ' 44 -0.22 


. .,-+0.20 
107 -0.29 


n „-+0.12 

°- 27 -o.io 


. o „+0.21 
1 -^'-0.09 


n - o +0.19 
07S -0.16 


. ,, o + 0.20 
1U8 -0.26 


Nh 




< 19.57 




4 40+ 1 - 57 




, co + 2.64 
3.»»_i. 96 


kTj 


, fi _+0.82 
^'"'-0.82 


n Rc+013 
U - DO -0.04 


J ' u2 -0.42 


46 +0 - 14 
u -™-0.14 


7 69+ 7 ' 08 
'■ Dy _3.05 


„ ^q + 0.05 
u - oa -0.07 


kT 2 




1 63+ 2 26 




1 - M1 -0.29 




4 27+ 4 ' 81 
^• z '-2.39 


T 


0.60+°;^ 




■| in + 0-51 
-0.35 




n O5+0-07 

u '^°-0.05 




EMi 


> 8.80 X 10 4 


55 +25 ' 35 
u ° -0.17 


03+ 62 
u - u,3 -o.oi 


R ,,+57.90 


< 0.44 


<270.41 


EM 2 




< 45.06 




o.82i° ;f 4 




<45.68 


He 


493 4R+ 9492 - 43 


< 202.72 


5319.53t™°| 5 94 


< 171.03 


„_ _+9580.50 
J ' a ' 1 ' -145.95 


PS 0.00 


C 


< 9999.98 


a 0.00 


a 0.00 


< 445.56 


pa 0.00 


a 0.00 


N 


SS 0.0 


PS 0.00 


< 9999.53 


< 805.06 


fa 0.00 


PS 0.00 


O 


4947. 23+l S ^ S1 


> 137.01 


< 100.83 


> 137.33 


< 421.57 


PS 0.00 


Ne 


< 195.20 


217.83t; 22 7 »« 


< 200.47 


o, -,+ 103.68 


< 45.27 


< 3.95 


Mg 




84 - 33 t37.95 




79.94+^"? 2 


64.24+115, 3 / 


„- , 7 +29.04 
^'■ 1 '-0.64 


Si 


256 53+ 52 64 
ZODM -52 .64 


101 23+ 54 24 
1UJ - M — 31.93 




163.73± 92 ' 9 7 


80.84™ 


6 40+ 29 71 


S 


298 4fi+ 105 - 28 


191 89 + 81 " 16 
j-»j-.o»_ 69 29 


1090.01+f 7 f;f 7 


125.29i« o ;™ 


1 93 7R + 277.87 


oe+207.92 
21.8o_ 1 64 


Ar 


SS 0.0 


< 73.88 


si 0.00 


101.34™° 


0.00 


> 1.32 


Ca 


> 862.33 


< 318.32 


< 3889.83 


< 262.68 


< 948.97 


< 5.50 


Fe 


780.96± 2 2 ^ 


90-78l S 4; 92 


< 237.873 


o fi 9 , + 79.64 


104 73+ 202 - 96 


- ,„ + 50.63 
°- d8 _1.07 


X?(d.o.f.) 


1.06(164) 


1.01 (162) 


1.44(173) 


1.32 (171) 


0.98(89) 


0.95 (87) 


Fx 


0.23 


0.23 


0.26 


0.26 


0.19 


0.19 


T T 




32.391 




32.451 




32.309 


L X 




31.941 




31.959 




31.687 


T H 




32.200 




32.282 




32.190 




Note: Hero Xti = X 2 


where v is degr 


ees of freedom. Errors are with 90% 


confidence based on 


X 2 mtn +2.71. 
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3.3 Evolution of spectral parameters 

The spectral analysis at different phases of the WR binaries 
provides the dependence of the best fit values of parameters 
on orbital phase. The variation of the soft (L x ) and hard 
(Li) band X -ray luminosities, column densities correspond- 
ing to the cool (N H ) and hot (Njj) temperature components, 
and cool(kTi) and hot (kT2) temperatures as a function of 
orbital phases of V444 Cyg and CD Cru are shown in Fig. 
[6ja) and[6jb), respectively. In the case of V444 Cyg, the 
Li was found to be minimum during the primary and sec- 
ondary eclipses and maximum outside the eclipse. It was 
constant during the phase 0.13 to 0.47. However, the L x 
was found to be minimum during the primary eclipse only. 
At the phase of 0.5 the L x was found to be maximum. Both 
Li and L x were found to be minimum during the phase 
0.0 and maximum at phase 0.5 in the WR binary CD Cru. 
The cool component (kTi) was found to be constant with a 
mean value of 0.61±0.05 keV and 0.57±0.10 keV for V444 
Cyg and CD Cru, respectively. The hard energy component 
(kT2) was varied from a minimum value of 1.88 keV to a 
maximum value of 9.61 keV for V444 Cyg. It was maximum 
at primary eclipse and at the phase of 0.29, and minimum 
during the secondary eclipse and at the phase of 0.13. For 
CD Cru, kT2 varies from a minimum value 1.63 keV dur- 
ing outside the eclipse to a maximum value of 4.27 keV at 
primary eclipse. For the star V444 Cyg, N H was maximum 
at phase 0.29, otherwise it was constant at all phases. How- 
ever, Njj was found to be constant throughout the orbital 
phase with a mean value of 1.14±0.11 xlO 22 cm -2 for CD 
Cru. For V444 Cyg, N| was increased from phase 0.0 to its 
maximum value at phase 0.13. Afterward it was decreased 
and became almost constant from phase 0.29 to phase 0.47. 
For CD Cru, Nh was minimum at phase 0.0 and maximum 
at phase 0.78. 



3.4 RGS spectra 

We combined the first and second order spectra of detector 
RGS1 and RGS2 to inspect the main spectral lines. We ob- 
tained the RGS fluxed spectra using task rgsfluxer and 
is shown in Fig. [7] for different phases for both the binaries. 
We identified the prominent emission lines for V444 Cyg at 
wavelengths Mg XII (6.50 A; 6.74 A), Mg XI (9.23 A), Fe 
XVIII (11.40 A) and Ca XVI (22.61 A, 21.43 A, 22.30 A). 
Individual lines show intensity variations from spectrum to 
spectrum. Unfortunately, our observations are at the limit of 
detectability with the data available with poor count statis- 
tics. Therefore, the signal-to-noise ratio of individual lines 
are not sufficient to perform a detailed quantitative analysis 
of the RGS line spectrum for any of the massive stars. 
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4 DISCUSSION 

We have analyzed the X-ray emission properties of two WR 
binaries with strong stellar winds. We discuss below the im- 
plications of the X-ray results. 



Figure 6. Variation of L|, Lg, Njj, N^, kTi and kT 2 as a func- 
tion of orbital phase of the the stars (a) V444 Cyg and (b) CD 
Cru. 



4.1 Phase-locked X-ray variability 

The X-ray temporal and spectral analyses of WR binaries 
V444 Cyg and CD Cru showed a phase-locked phenomenon. 
The phase resolved X-ray spectroscopic observations show 
that the soft energy component and the corresponding N H 
were also found to be nearly constant throughout the binary 
phases of both stars. The L x was found to be minimum dur- 
ing the primary eclipse only. These results lead to the ra- 
diative wind shock origin of soft energy component. In the 
case of CD Cru, the modulation of L x well matches with the 
optical light curve reported by Moffat et al. (1990). Lamon- 
tagne et al. (1996) have reported a similar shape of optical 
photometric light curves for 13 WR binaries including CD 
Cru and they said that it could be due to the "atmospheric 
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Figure 7. RGS spectra of the WR binaries at different binary phase for the stars V444 Cyg (upper panel) and CD Cru (lower panel). 



eclipse" . The X-ray flux in the hard energy band was found 
to be minimum at both primary and secondary eclipse for 
both stars. This implies that the hard energy component 
could be due to the presence of wind-wind collision zone 
which is located somewhere in between the O-type and WR 
star, and therefore, minimum when either of the star in the 
binary system is eclipsed. A similar behavior of X-ray light 
curves was seen in the EINSTEIN, ROSAT and ASCA ob- 
servations of V444 Cyg (Moffat et al. 1982; Corcoran et al. 
1996 ; Maeda et al. 1999). The hot temperature (kT 2 ) and 
the corresponding Ng were also varied with orbital phases 
of both binaries. The eccentricity of both binary systems 
is almost zero, therefore, the variation in the temperature 
corresponding to the hard energy component could be the 
result of variation in N^, i.e., varying optical depth (see SJTJ). 
This further supports the wind-wind collision phenomenon 
in these binaries. 



4.2 X-ray Temperatures of Plasma 

The X-ray spectra of WR binaries are well defined by two 
temperature plasmas. The values of log(Lx/Lboi) in the total 
energy band are found to be -6.6 and -7.2 in for V444 Cyg 
and CD Cru, respectively, which are similar to those for 
other WR (WN) stars observed from XMM-Newton and 



CHANDRA (Gudel & Naze 2009). The temperature of the 
cool component of WR binaries are comparable to other 
WR binaries e.g., 0.56-0.67 keV for WR1 (Ignace, Oskinova 
& Brown 2003) and 0.7-0.8 keV for WR 147 (Skinner et 
al. 2007). However, the temperatures corresponding to hot 
energy component o f these WR binaries are slightly more 
than some of the similar binaries. The possible mechanisms 
for generation of X-rays are discussed below. 



4-2.1 Instabilities driven radiative wind shocks 

The wind-shock model predicts the intrinsic instability of 
the line driving mechanism. The standard model estimates 
the shock velocities by the relation kT s h = l-95^tVg hock (Lucy 
1982; Luo, McCray & Maclow 1990). The temperatures of 
the cool energy components are found to be almost simi- 
lar (~ 0.6 keV) for both the binary systems (see Table [3] 
and Table [4]). Adopting the mean particle weight n ~1.16 
(Skinner et al. 2007) for WN stars and fx Ri0.62 for O-type 
stars (Cassinelli et al. 2008), we derived the "average" value 
of shock velocities « 515 km s _1 for WN stars and « 704 
km s _1 for O-type stars corresponding to the cool compo- 
nent. These values are about a factor of 2 larger than those 
predicted by radiative shock model of Lucy (1982). However, 
the advance version of the wind shock model by Owocki et 
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al. (1988) predicts X-ray emission up to 1 keV. Moreover, 
the temperatures of the cool energy component and Njj are 
not varying with orbital phase. Therefore, it appears that 
the cool energy component from V444 Cyg and CD Cru 
could be generated by either of the binary components and 
may be explained by instabilities in radiation-driven wind 
shocks. 



4-2.2 Magnetically Confined Wind Shock 

Babel & Montmerle (1997) have suggested that the presence 
of magnetic fields confine the wind, which may be an impor- 
tant ingredient for the production of X-ray emission (Babel 
& Montmerle 1997). The degree of confinement of wind by 
the magnetic field is derived in terms of a confinement pa- 
rameter r = BqR^/Mvoo, where Bo is the surface equatorial 
magnetic field strength (ud-Doula & Owocki 2002). For a 
confined wind model T >> 1. Using the values of stellar ra- 
dius (R») and stellar mass loss rate (Mv) as given in Table [1] 
for V444 Cyg binary components, the minimum magnetic 
field required to confine the wind is derived to be 0.14 kG 
for the primary and 1.25 kG for the secondary. Similarly, 
for CD Cru, we derived minimum magnetic field of 0.16 kG 
and 1.96 kG for the primary and the secondary components, 
respectively. Based on the formula given by Babel & Mont- 
merle (1997) and using above estimated values of B, the 
X-ray luminosity are estimated to be io 34 - 66 erg s _1 and 
10 35 ' 74 erg s^ 1 for primary and the secondary components 
of V444 Cyg, and 10 35 06 erg s" 1 and 10 36 8 erg s" 1 for 
primary and the secondary components of CD Cru. Such 
a high X-ray luminosity has not been observed from any 
of the WR binaries. However, smaller fields can exist and 
the confinement may be limited to just above the magnetic 
equatorial plane. 



4-2.3 Colliding wind shock model 

Assuming that shock velocities have reached to the terminal 
velocities, the standard model predicts a maximum temper- 
ature for the shocked region from winds of massive stars 
by the relation kT^ ax = 1.95/w^ (Luo, McCray & Maclow 
1990). The nominal wind parameters of massive stars (see 
Table [1]) and the mean particle weight fj, for the WN stars 
and for the O-type stars give the maximum temperature gen- 
erated by the individual components of binaries. The derived 
values are 7.80 keV + 7.21 keV for V444 Cyg (06+WN5) 
and 10.88 keV + 13.69 keV for CD Cru (05+WN6). The 
maximum observed temperature at phase 0.29 correspond- 
ing to hot component is found to be similar to that of max- 
imum possible values in the case of V444 Cyg. However, in 
CD Cru maximum observed temperature is lesser than that 
of maximum possible values. 

The phase locked variability of the hard energy com- 
ponent shows that the observed hot X-ray temperatures in 
massive binaries could originate from the collision of stellar 
winds. The distances from the stars where these winds meet 
are derived using the relations (De Becker 2007; Stevens, 
Blondin & Pollock 1992), 



where tob and D are the distance to the collision zone from 
the primary and the separation between the binary compo- 
nents, respectively. The wind momentum ratio r\ is expressed 
as, 



V = 



M 2 v D 



Miv 



(2) 



where Mi and M2 are mass loss rates of primary and sec- 
ondary components, respectively, and Voo,i and v-x,^ are ter- 
minal velocities of primary and secondary components, re- 
spectively. Using the above relations, we derived wind mo- 
mentum ratios of 7.03 and 24.60 for the stars V444 Cyg 
and CD Cru , respectively. Distances to the collision zone 
from the primary are estimated to be 0.27 and 0.17 times 
of their corresponding binary separation (D) for V444 Cyg 
and CD Cru, respectively. Therefore, the collision zone ex- 
ists to be very close to the primary O-type star in case of 
both the binaries. Further, the collision wind zone will pro- 
duce a bow shock around the O-type star , the shock volume 
and the emission measure are dominated by the WR winds, 
therefore, X-ray emitting plasma also shows non-solar abun- 
dances. 

The gas in a colliding wind region could be either adi- 
abatic or radiative and depends upon the value of cooling 
parameter (\), which is defined as the ratio of the cooling 
time (t coo i) of the shocked gas to the escape time (t osc ) from 
the intershock region (Stevens, Blondin & Pollock 1992). 



X = 



v 3 d7 

M-7 



(3) 



where V3 is the pre-shock velocity in the unit of 10 3 km s , 
d7 is the distance from stellar center to the shock in the 
unit of 10 7 km and M_7 is the mass loss rate in the unit 
of 10 _7 MQyr _1 . Using the parameters given in Table[T]and 
the values from above estimation, we derived cooling pa- 
rameters xi an d X2 f° r primary and secondary components, 
respectively. The estimated values of xi an d X2 are given 
Table The winds from WR stars are found to be clearly 
radiative (x < 1) for both the binaries. The intrinsic lumi- 
nosity can be estimated for each component of the binary 
using the following relation (Pittard & Stevens 2002) : 



L x = 0.5 3 M v 



(4) 



roB = 



,1/2 



D 



(1) 



where 3 accounts for a geometrical and inefficiency factors 
and v is the wind speed at contact surface (i.e. pre-shock ve- 
locities derived from observed temperatures). The values of 
Hi (0.403 for V444 Cyg ; 0.564 for CD Cru) and H 2 (0.033 for 
V444 Cyg ; 0.0042 for CD Cru) for primary and secondary 
stars, respectively, are taken from Pittard & Stevens (2002). 
The estimated values of intrinsic luminosities of each com- 
ponent and the total intrinsic luminosities of binary systems 
are given in Table [5] These values are nearly 3 orders of mag- 
nitude higher than those of the observed values. A similar re- 
sults have been found for other WR binaries, e.g., HD 159176 
(De Becker et al. 2004); WR 147 (Skinner et al. 2007). De 
Becker et al. (2004) suggested that the disagreement be- 
tween observed and theoretical predictions could possibly 
be explained by (a) the kinetic power of the collision should 
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Table 5. Derived values of wind parameters and intrinsic luminosities for binary systems. The pre-shock velocities for each components 
are estimated using derived values of post-shock temperatures. 



Object 
Name 



Phase 



* pre — shock 

km s — 1 



log( L X i) 
,„._ „— i 



log(L X 2) 
„.„ „-i 



log(L total ) 
erg s _1 



V444 Cyg 



0.01-0.04 
0.12-0.14 
0.27-0.30 
0.38-0.41 
0.45-0.50 
0.50-0.51 



2598 

1237 
2819 
1735 
1632 
1276 



745 
779 " 
+71 
90 
475 
306 " 
+226 

249 " 
+164 

89 
+ 1242 
358 



H 904 ±66 

■ 1269^1 

- 1193tg5° 
+ 033+9O8 



5.496 


0.416 


35.711 


0.282 


0.021 


35.067 


7.623 


0.577 


35.782 


1.094 


0.083 


35.361 


0.856 


0.065 


35.307 


0.320 


0.024 


35.094 



+0.109 
-0.155 
+0.024 
-0.033 
+0.068 
-0.050 
+ 0.053 

-0.067 
+0.042 
-0.024 
+0.295 
-0.143 



35.352 
34.708 
35.423 
35.002 
34.948 
34.735 



+0.109 
0.155 
+0.024 
0.033 
+0.068 
0.050 
+ 0.053 
0.067 
+0.042 
0.024 
+0.295 
0.143 



35.869 
35.224 
35.940 
35.518 
35.465 
35.252 



+0.087 
0.109 
+0.021 
0.022 
+0.043 
0.048 
+0.043 
0.048 
+0.025 
0.026 
+0.169 
0.281 



CD Cru 



0.01-0.09 
0.43-0.51 
0.73-0.83 



1879 
1224 
1161 



+861 
632 
+ 101 
-102 
+633 
85 



1374±«0 

f 896±?S 
-849±f 2 2 



0.990 
0.178 
0.144 



0.047 
0.008 
0.007 



35.7981°;$ 



35.425~* 



-+0.034 
O'lll 

35.380_ 033 



34.8751°;!?! 
34.502t°;°§t 
34.4571^1 



35.847lg;ig 
ao.iii_ 034 



35.428 



+0.107 



35.8 



35.6 



m 35.4 

M 

<U 



36.0 
35.8 
35.6 
35.4 
35.2 
35.0 



i 1 

CD Cru 



V444 Cyg 



32.1 32.2 32.3 32.4 32.5 32.6 32.7 

log(LJ) (erg s" 1 ) 

Figure 8. Observed vs Model X-ray luminosity. 



be considered as an upper limit on the X-ray luminosity, (b) 
higher value of the parameter r\ should be considered, (c) 
diffusive mixing between hot and cool material is likely to 
exist due to the instability of the shock front, and (d) orbital 
effects should also be included to study such systems. The 
observed and theoretically predicted values of X-ray lumi- 
nosities are shown in Fig. [8] This model also predicts the 
phase-locked variations which can be seen in Fig. [6] There- 
fore, it appears that the hard energy component is most 
likely associated with wind collision zone. 



4-2-4 Non-thermal emission 

There is no strong justification for invoking non-thermal ori- 
gin of X-ray emission since the observed X-ray emission at 
higher energies can be fitted satisfactorily with a thermal 
model (see Fig. [4] and Fig. [5}. However, the present data 
is limited up to 10 keV and it is difficult to eliminate the 



contribution of non-thermal components from the spectra 
which is dominated by the presence of the strong thermal 
emission, and the non-thermal emission is overwhelmed by 
the thermal emission (De Becker 2007). Therefore, we can 
state that the spectra of massive stars in the present sample 
are fairly well explained by two-temperature plasma models. 



5 SUMMARY 

We have analysed X-ray temporal and spectral properties of 
two WR binaries V444 Cyg and CD Cru observed with high 
sensitivity EPIC instruments on board the XMM-Newton 
satellite. The X-ray light curves in the soft and the hard 
energy band show the phase locked variability. Both pri- 
mary and secondary minima were seen in the hard band 
X-ray light curves of both binaries. However, in the soft 
X-ray light curve only primary minimum was seen. This im- 
plies that the hard energy component could be originating 
from the wind- wind collision zone. The X-ray spectra of both 
WR binaries show strong absorption below « 1 keV and a 
clear evidence of high temperature plasma, manifested by 
a visible Fe Kct emission-line complex. The study of X-ray 
spectra reveals cool as well as hot temperature plasma com- 
ponents of binary stars which are fitted consistently with 
two-temperature plasma models. The cooler plasma compo- 
nent was found to be constant at all phases with mean value 
of ~0.6 keV for both binaries. The presence of cooler compo- 
nent could be due to the distribution of small-scale shocks in 
the radiation-driven outflows from either the primary or the 
secondary star in the binary systems. The temperature of 
hot plasma component and the corresponding column den- 
sity were found to be variable during the orbital cycle of both 
binaries. The variation in temperature of hot plasma could 
be due to the varying circumstellar optical depth along the 
line of sight towards the shock as star revolves around each 
other. The maximum value of hot plasma was found to be 
lower than the hottest plasma possible in the binary systems 
as predicted by colliding wind theory for short-period bina- 
ries. However, the predicted values of X-ray luminosities are 
~ 3 orders of magnitude more than those of the observed 
values and can not be accounted in terms of observational 
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